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Abstract 

The approach to equilibrium for systems interacting with their environment 

by being regularly exposed to low energy, low intensity pulses of some type 

of quanta is studied. Assuming a randomness condition on the interaction of 

these quanta with the system, but making no assumptions about the accessible 

states of the system, we show that microcanonical equilibrium is reached. 

Although the intensity of the pulses is assumed to be weak the interactions 

are treated exactly, with no recourse to perturbation theory. The origin of 

time asymmetry and the absence of recurrence is discussed. 
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I. INTRODUCTION 



For over a 100 years various approaches have been used to study the foundations of 
statistical mechanics. That this subject continues to be of interest |lj indicates that no 
completely satisfactory description exists; this is true, both, for classical and quantum sta- 
tistical mechanics. The question of how macroscopic irreversibility arises from reversible 
microscopic dynamics continues to be of a topic of discussion and is reviewed in Ref. ffl. 

For quantum statistical mechanics the goal to show that a system can be described by a 
microcanonical density matrix, which for some given energy E is 

where N(E) is the number of degenerate eigenstates with energy E. It would be desir- 
able to obtain the completely mixed state represented by the density matrix of eq. ([[]) as 
a time development of a pure initial state; this is of course impossible under unitary time 
evolution and one has to resort to various averaging assumptions. In the most straight for- 
ward ensemble "derivation" of this result one essentially assumes the answer by postulating 
random phases and equal d priori probabilities for the eigenstates of some Hamiltonian. 
More recently developments in the theory of chaotic systems have lead to a different ways 
of obtaining quantum statistical mechanics [fj . An old approach to ([I]) is through the Pauli 
master equation ||, 

= E [WnW " tn^iVm , (2) 

where Vi(t) is the probability at time t of the system being in the state i and is the rate 
per unit time to go from state n to state i. In order to obtain eq. (§) certain assumptions 
have to hold [[|: 

(i) A repeated randomness assumptions insuring that the off diagonal elements of p { j 
vanish. This is an assumption on the states that are accessible as the system evolves. 
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(ii) The interaction potential that causes the mixing of states is assumed to be weak as to 
permit the use of first order perturbation theory. 

(iii) There is an inherent clash in that the TVs refer to discrete states and yet one has to 
use continuum normalization in order to obtain a Dirac delta function in the energies. 

Conditions (i) and (ii) can be weakened for certain classes of interactions Q (also see L. Van 
Hove's in Ref. [|]). 

In this presentation we shall consider a nonisolated, finite but large, system that is 
allowed to interact with its environment; by this we mean that from time to time the system 
is exposed to a low energy, low intensity packet of some quanta. These quanta are spatially 
unconstrained (continuum normalized) and after a certain characteristic interaction time 
r the system settles into a new state or density matrix. After intervals larger than this 
characteristic time, this process is repeated. The above qualitative terms, low intensity, low 
energy and time interval, will be made precise and summarized in Sec. |V[ The quanta in 
these pulses will be referred to as "photons". We shall show that with one assumption on 
the interaction of these packets with the system, but independent of any assumptions on 
the accessibility of the system itself, microcanonical equilibrium, eq. ([!]), will be reached. 
Contrary to the approach to equilibrium for isolated systems, in this case equilibrium is 
attained by "entangling" the system with photons and then looking at expectation values 
of operators that are diagonal in the photon variables. Even though we may start out with 
a pure system-photon state, the system itself will be in a mixed state after the interactions 
have ceased. No recourse is made to perturbation theory; the interactions are treated exactly. 

The system, its interactions and time evolution are presented in Sec. [II] while the time 
evolution of density matrices is discussed in Sec. [TIT]; an assumption on the randomness of 
phases of certain amplitudes is presented and discussed in this section. How the system 
approaches equilibrium is shown in [TV]; although we do not get a master equation the 
approach to equilibrium is Markovian [H] , governed by a Chapman-Kolmogorov [|3| equation. 
The relations among various energy and time scales that the system and perturbing quanta 
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have to satisfy are given in Sec. [V| where the origin of time irreversibility and the absence 
of recurrence is discussed. 



II. INTERACTIONS AND TIME EVOLUTION OF STATES 



The energy levels E a of the system are highly degenerate with states \a, i) .The energies 
of the photons in the packets are taken to be too small to cause transitions between states of 
different energy; these interactions can, however, cause transitions between the degenerate 
levels within a given E a \ for convenience we take the energy E a = and drop the quantum 
number a in the description of the states. As usual the total Hamiltonian will be split into 
two parts, H = Hq + Hi; eigenstes of Hq that are of interest are \i, k) with 



the specific dispersion relation for uj(k) is unimportant. Hi induces transitions between 
\i, k) <-» \ j, k'). Whether this Hamiltonian is time reversal invariant or not is immaterial to 
subsequent developments. 

Let \i,k)fj be that eigenstate of the total Hamiltonian that approaches \i,k) for large 
negative times; we use the notation | )# for an eigenstate of the total Hamiltonian, whereas 
| ), without the subscript H, for an eigenstate of H Q . This state satisfies the Lippmann- 
Schwinger equation || 



H \i, k) = u(k)\i, k) ; 



(3) 



i, k) H = \i, k) + 



1 



H 1 \i,k) H . 



(4) 



u{k) — H + ie 



The overlap with an eigenstate of the free Hamiltonian is 



(j,p\i,k) H = S itj 8(k-p) + 



1 



(5) 



ui{k) — ui{p) + ie 



where 




(6) 
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is the scattering amplitude for the transition \i, k) — > \j, p)and satisfies the off shell unitarity 
relations 

1 f 



T J>l (p,k)-T* ] (k,p) = Jdq 

dq 



oo(q) — u(p) + ie u(q) — u(k) — ie 
1 1 



T jin (p,q)7^ n (k,q), 

T nj (q-k)T; Jq.p) . (7) 
uj{q) — uj(p) + ie uj{q) — u(k) — ie\ ,J 

We ant to describe the time evolution of a state that at t — started out as |i, fc); more 
specifically we want the overlap at time t with the state \j, p). 

(j, p\e- lHt \i,k) = [5 id 6(p- k)+iTZ j , i (p,k;t))e- i ^ t ] 



(8) 



in the above 

iTZj,i(p,k;t) 



u{p) — tu(k) — ie 



-r^p, ky^-^ + r^p) 



+ / dqe 



i[u)(p)~u>(q)]t _ 



Tj,n{p,q) 



(9) 

u(q) — u(p) + ie u(q) — u(k) — ie 
For the definition of 1Z in eq. (||) we have pulled out the factor exp [— iuj(p)} for later conve- 
nience. Eq. (|7]) or (||) implies unitarity relations for 1Z 

- i [K jti {k', k; t) - ni 3 (k, k'; *)] = E / dqn j>n (k', q, tR* n ( fc, q, t) , 

n J 

= E / dqK n:i (q, fc; t)H* n /q, k'; t) . (10) 



III. TIME EVOLUTION OF THE DENSITY MATRIX 

We shall be interested in operators O be that have off diagonal matrix elements between 
the different |i)'s and are diagonal in the photon subspace. For any state 

\S) = J2jdqMq)\hq), (11) 

i 

the expectation value of O is 

(S\0\S) = 52(j\0\i) Pi>j , (12) 

with a density matrix 

p w = UqUqWM)- (is) 



5 



A. Evolution of p entangled with a wave packet of photons 



Suppose that at time t = our system, described by a density matrix p -^O), is exposed 
to a photon state. The total density matrix, for the system plus photons is 

P T = Y,j dkdk^(k)r(k')\j,k) Pj>i (0)(hk'\ , (14) 

where ip(k) describes the photon wave packet scattered of the mixed state at t — 0; 
/ dk\%l){k)\ 2 = 1. Using eq. and eq. (^) the density matrix at time t, after summing 
over the photon states, is 



p Jt i(t) = Pj,M + / *^(*)f(*')[iE^(*'.*;*)Pn 1 i(o) 

J n 

-iJ2Kl m (k,k>;t)p jm (0)+ f dp^7e i>n (p,fe;t)^* ni (p ) fe / ;t)p„ )m (0) 

m n,m 

(15) 

In subsequent discussions we shall need some randomness condition on the phases of the 
72.'s. It is unlikely that such a condition could be valid for all times; we shall try for ones that 
may hold at large times. For short duration pulses we expect that after some characteristic 
collision time r, as for example the inverse of the width of a resonance in resonance dominated 
scattering 0, the system will settle down and 

lim f dki){k)Kjj{p, k, t) = K hi (p) ; (16) 

lZj t i(p) depends implicitly on the wave function ip and, using the definition in eq. @ we 
find 

n iAP ) = jimk )w ^±^— Wp) . (17) 

For further developments, this explicit form will not be needed. We also define 

n jti = J dtep*(k')K jti (k') . (18) 
These, in turn, satisfy the unitarity relations 



- i (n hi - 72*,.) = £ / dpKnMKM • ( 19 ) 

In particular, we find 

Im72 iJ >0. (20) 
The evolution of the density matrix, eq. (|i~5l) , may be expressed as 

Pj,M = + 4 £ [^,«Pn,<(0) " ^nPj>(0)] + £ / dP^,«(p)^(p)Pn,m(0) • (21) 

n m,n 

B. Randomness Assumption 

In order to proceed further we must impose a crucial condition on the 72's: 

J dpK j)n {p)K\ t Jp) = for i ^ j and m ^ n. (22) 

This results from the assumption that, as we integrate over p, the phases of TZj tn (p) fluctuate 
rapidly. It should be emphasized that this is an assumptions on the dynamics of the system 
and not on states or density matrices of the system at any particular time. 
The above and unitarity relation, eq. (0), lead to 

(ll hl - 111) =0 for i^j; (23) 

which together with eq. Q2Tj) implies 

K = K H + iV (24) 

with TZh Hermitian and T> a diagonal matrix with positive elements. As a matter of fact, 
there always exists a basis of the states \i) where such a decomposition of 72 holds. Any 
matrix can be written as a sum of a Hermitian and an anti-Hermitian one and we go to 
the basis where the anti-Hermitian part is diagonal. Although eq. (|22|) implies eq. (0), the 
inverse is not true and eq. (^) remains an assumption. We work in the basis where this 
decomposition holds. The time evolution of the density matrix becomes 

Phl (t) = Pjti (o) + 1 [n H , P (o)]^ - [v, p(o)] u + Stj £ / d P \n lin (p)\ 2 PnA°) ; ( 25 ) 

[A, B] is the commutator of the matrices A and B. 



IV. APPROACH TO EQUILIBRIUM 

At this point we have to require the IZj/s to be small; how this is achieved will be made 
precise in Sec. M. We solve eq. (P5f) to first order in the IZ's by first rewriting it as 

Pj,M = Pj,M - P> p(°\i + E / d P \n ltn (p)\ 2 PnA°) > ( 26 ) 

n J 

with 

hA 1 ) = E [i - p n>m (t) [i + in) m . . (27) 

n,m 

To the order we are working 1 — iTZ is a unitary matrix and Pji(t) is the density matrix in 
a basis rotated from the one we started out at t — 0. 
Let us first look at eq. fl26| ) for i ^ j. 



p h ,(l) 1 CP,.; • !>,,((>). (28) 

The off diagonal elements of the unitarily rotated density matrix at time t are smaller than 
the corresponding matrix element at t = 0; 

For the case i = j we have 

Pm(*) = + * fa* - n t) M°) + E / dp\Ki, n (p)\ 2 p n>n {0) ■ (30) 

The coefficient of p nn (0) in the last term of this equation may be identified with the tran- 
sition probability for \n) to evolve into 

W^ n = [ dp\U hn {p)\ 2 . (31) 



Using the unitarity relation, eq. (|19|), we find 

i (n hl -ni l ) = -J2^n^ (32) 

n 

and the evolution equation for this case becomes 



P M (t) = p ifi (0) - £ W n ^p hl (0) + £ wu^(o) • (33) 

n n 

Even if we include the unitary rotations this equation is of the Chapman-Kolmogorov type 
@ and for the Wj^„'s not too large, will directly yield microcanonical equilibrium. This 
follows from the observation that the matrix 

M iJ = -6 i>j J EWn~i + Wij (34) 

n 

has one eigenvalue equal to zero and all others are negative which in turn is obtained by 
showing the function H{r) = — J2i'Pi{ T ) lnT^r) satisfies a Boltzmann H-theorem, dH/dr > 
with the probabilities being functions of the auxiliary variable r and (c.f. eq @) 

f = EM/ r (35) 

Let ~ (1, 1, ■ • • , 1) be the eigenvector of with eigenvalue zero and be all the 
others. 

Pi,i = CQvf ] + £ c a v\ a) . (36) 

a 

At the end of the interval t cq hasn't changed and the other c Q 's have all decreased in 
magnitude. We find that repeated pulses will drive off diagonal elements to zero and the 
diagonal ones to the same constant or the final density matrix will be as in eq. ([[]). 



V. CONCLUDING REMARKS 

We have to consider three characteristic energies: (i) A a , the difference in the energy 
levels E a of the system, (ii) 5 W , the energy spread of the impinging wave packets, and (iii) 
7, the inverse of the interaction time r. These quantities have to satisfy 

A a > 5 W > 7 . (37) 

The first of these inequalities insures that there will be no transitions between the different 
energy levels of the system; for such transition the energy denominator in eq. @ would have 
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a typical magnitude of A a as opposed to 5 W for intra level transitions. Due to the second 
inequality the "arrival period" of the pulse is much shorter than the interaction time and 
for times greater than r the system and the photons propagate separately. In addition the 
repetition time of these pulses, T, must satisfy 

T > - ; (38) 

7 

this guarantees that as a new packet interacts with the system it is unencumbered by photons 
from the previous pulses; for example T must be large enough to allow any possible system- 
photon resonances to decay. 

Till now, all this has been worked out for the system interacting with a one photon packet; 
the random impulses from the outside are likely to contain many photons and we extend 
our results to iV photons where the states are / IIn=i <^k n ip n (k n )\i; ki, k 2 • • • k/v). In the 
case the packets ^n(k) are different for different n's this expression should be appropriately 
symmetrized; doing this explicitly would just lead to unnecessary notational complications. 
The scattering amplitude is not necessarily a sum of two body amplitudes and the rest of 
the development is as the previous one. We identify a weak, or low intensity pulse as one 
with few photons and an intense one with a large number. We require that N be sufficiently 
small that the TZj/s are small compared to one; N should not be too small as the time to 
reach equilibrium would become very large. 

With these conditions and the assumption about interactions made in eq. ( j22|) satisfied, 
repeats interactions of the system with external quanta will drive it to equilibrium. Time 
irreversibility is built in right at the start in the choice of the sign of the ie term in eq. (f|) 
The continuous energy levels of the photons preclude any recurrences; had the photons also 
been quantized in a finite volume the limit considered in eq. (^6j) would not have existed. 
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